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We elucidate the geometry of quantum adiabatic evolution. By minimizing the deviation from adiabaticity 
we find a Riemannian metric tensor underlying adiabatic evolution. Equipped with this tensor, we identify a 
unified geometric description of quantum adiabatic evolution and quantum phase transitions, which generalizes 
previous treatments to allow for degeneracy. The same structure is relevant for applications in quantum infor- 
mation processing, including adiabatic and holonomic quantum computing, where geodesies over the manifold 
of control parameters correspond to paths which minimize errors. We illustrate this geometric structure with 
examples, for which we explicitly find adiabatic geodesies. By solving the geodesic equations in the vicinity of 
a quantum critical point, we identify universal characteristics of optimal adiabatic passage through a quantum 
phase transition. In particular, we show that in the vicinity of a critical point describing a second order quantum 
phase transition, the geodesic exhibits power-law scaling with an exponent given by twice the inverse of the 
product of the spatial and scaling dimensions. 
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I. INTRODUCTION 

Geometric and topological concepts have long played use- 
ful roles in both classical and quantum physics [IJ. Important 
applications where the use of geometry has led to new insights 
include quantum evolutions |2|, distance measures in quan- 
tum information theory |3 , 4|, circuit-based quantum compu- 
tation Is], and holonomic quantum computation ||6|. More 
recently quantum phase transitions (QPTs) |7| and adiabatic 
quantum computation |8 9| have also been been explored 
from a geometric perspective iflOl fTTl . While geometry can 
be seen as an underlying unifying theme in these applications, 
an explicit geometry-based connection between them is not 
always apparent. The central theme of this work is to elu- 
cidate the geometry of adiabatic evolution. In particular, we 
describe an all-geometric connection between QPTs and adi- 
abatic quantum evolution. We do this by showing how the 
Riemannian metric tensor that describes transitions through 
quantum critical points IITOl also arises in adiabatic quantum 
evolution. More specifically, we explain how the metric which 
provides an information-geometric framework for QPTs can 
also provide a geometry for the control manifold arising in 
adiabatic evolutions. That QPTs and adiabatic quantum evo- 
lution should be so intimately related was previously under- 
stood in terms of the role of ground state evolution in adiabatic 
quantum computation, and in particular the basic observation 
that those points where ground state properties undergo dras- 
tic changes, i.e., quantum critical points, are bottlenecks for 
adiabaticity ISl [Til IT3l . 

The metric tensor we identify is a natural extension of the 
metric found in Ref. fTU\ to systems with degenerate ground 
states. In this sense we go beyond adiabatic quantum com- 
putation, which is typically concerned with nondegenerate 
ground states, and find results with applications to holonomic 
quantum computation, where quantum gates are performed 
as holonomies in the degenerate ground eigensubspace of the 
system Hamiltonian. We analyze the relevance of the metric 



tensor we identify for determining paths with minimum com- 
putational error, in the sense of deviation from the desired 
final adiabatic state. In addition, we find a prescription for 
adiabatic passage through quantum critical regions by solving 
the coiTesponding geodesic equations derived from the metric 
tensor. As a result we are able to identify universal charac- 
teristics of adiabatic passage through a critical point. Namely, 
we find that in the vicinity of a critical point the geodesic ex- 
hibits power-law scaling with an exponent given by twice the 
inverse of the product of the spatial and scaUng dimensions. 

The structure of this paper is as follows. In Sec.|ll]we for- 
mulate our geometric picture. Specifically, after defining the 
model in subsection |II A[ in subsection |IIB| we introduce the 
adiabatic error and show how to upper bound it as a sum of 
two components, one of which encodes the geometric aspects 
of the evolution. We obtain a Riemannian metric by mini- 
mizing this error. Next, in subsection II C we demonstrate 



the emergence of the same geometry from the concept of adi- 
abatic operator fidelity. In subsection II D we demonstrate 



how our metric arises from three more (interrelated) natural 
origins: Grassmannian geometry, Uhlmann parallel transport, 
and the Bures metric. In subsection II E we compare our 



metric with another, related metric for adiabatic evolutions 
which we proposed in earlier work [11 1. We briefly discuss 
strategies for further making the adiabatic error small in sub- 
section InFl We make the connection to QPTs in section [lll| 



Specifically, in subsection III A we establish the relevance of 
our metric in the sense of QPTs, by showing that the same 
metric is responsible for signaling quantum criticality. Then, 



in subsection III B we derive the quantum critical scaling of 
the metric tensor Switching gears, we define the notion of 
an adiabatic geodesic in Sec. |IV] In subsection |IV A| we an- 
alyze three examples, namely the Deutsch-Jozsa algorithm, 
projective Hamiltonians (including Grover's algorithm), and 
the transverse field Ising model, for which we analytically find 
the adiabatic metric and the coiTesponding geodesies. In sub- 



section IV B we analyze the properties of geodesies when the 
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adiabatic evolution passes through a quantum critical point. It 
is here that we identify the universal characteristics of such 
geodesies. We summarize our results and conclude in SecfV] 
Several appendices provide detailed proofs omitted from the 
main text so as not to interrupt the presentation. 

11. GEOMETRY OF ADIABATIC QUANTUM EVOLUTION 
A. Model 

Consider an n-body system with the iV-dimensional Hilbert 
space H. The Hamiltonian family {_ff(x)} for this system, 
which depends on the (time-dependent) coupling strengths or 
"control knobs" x, can be identified by points over the real ill- 
dimensional manifold M 3 x. Given a total evolution time 
T and rescaled time s = t/T, a path x : s e [0, 1] i-> 
then represents the dynamics in this time interval, starting 
from Xq = x(0) and ending at Xi = x(l). We shall use 
the notation Xg = x(s) interchangeably, or sometimes drop 
the s-dependence entirely to lighten the notation. We al- 
low for a (7o(x)-fold degenerate ground-state eigensubspace 
of {H{x)}, with eigenstates {|<I>q (x))}. Thus this subspace 
can be identified by the projector 

90 (x) 

Po(x) = J2 1*0 W)(*o Wl. (1) 

Q = l 

with Tr[Po(x)] = 5o(x) > 1. We assume that for all finite n 
the ground-state energy Eo{x) is separated by a nonvanishing 
gap A(x) from the rest of the spectrum. In the thermody- 
namic limit n ^ cx) we allow the gap to vanish at some finite 
set of points {xc = x(sc)}, or a segment of the path. These 
are the critical points where a QPT takes place. Although our 
results would hold if we picked any other eigensubspace satis- 
fying the previous requirements, rather than the ground state, 
for specificity we shall henceforth consider the ground state 

and the initialization |V^(0)) = T,alT^ a^l'^oi^o)) (where 
\ip{s)) = \ip (x^)), and we similarly drop the explicit depen- 
dence on x(s) hereafter where possible). 

B. Adiabatic Error 

1. Degenerate case 

We wish to compare the desired, "ideal" adiabatic evolution 
to the actual evolution induced by the the Hamiltonian family. 
To this end we shall define an appropriate "adiabatic error" 
which measures the deviation between the two. The state of 
the system, 

ms)}=V{s)\m): (2) 

at any rescaled time s, is given in = 1 units [adopted here- 
after], in terms of the propagator V{s) which is the solution 
to the time-dependent Schrodinger equation 

idsVis) ^THis)Vis). (3) 



We can similarly associate an adiabatic propagator Kd(s) and 
an adiabatic Hamiltonian i?ad(s) to the ideal adiabatic evolu- 
tion, where the two are related via the Schrodinger equation 

idsV,a{s) = TH,a{s)V,a{s). (4) 

What defines the adiabatic propagator is the "intertwining 
property" 

Kd(s)Po(0)K|j(s) = Po(s), (5) 

which means that Kid(s) preserves the band structure of 
the ground eigensubspace of H{s). By differentiation 
the intertwining property is equivalent to idsPo{s) = 
T[Hai{s), Po{s)], and when it holds we have 

So 

\i^.,,{s))^V,,{s)\m)^ E «o^i°i'WI*o'(s)), (6) 

aa'—O 

where Fi^,(s) = ($^(s)|Kd(s)|$^'(0)) is the (non-Abelian) 
Wilczek-Zee holonomy (14] — usually expressed as the path- 
ordered exponential 

y[°l(s) =Pexp(-^ ^(s')ds'), (7) 

with the gauge connection 

A^^, = (8) 

We prove Eq. (|7| in Appendix[A|(see also Ref. ifTSi ). 

The adiabatic Hamiltonian can be expressed in terms of the 
original Hamiltonian plus a "correction" term lfT6l[T7l : 

ffad(s) - H{s) + i[dsPo{s),Pois)]/T, (9) 

Clearly then, the actual state IV'(s)) need not be the same as 
the adiabatic state |^ad(s))- Our objective is to find the path 
Xs that minimizes the adiabatic error || |'0(xs)) — |^ad(xs)) || — 
\\{V{xs) — T4d(xs)}|'!/'(xo))||, where the norm is the standard 
Euclidean norm: \\\(f))\\ = y/ ((/>|0). However, so as to obtain 
a result which does not depend on the initial state |?/'(xo)) we 
shall adopt a state-independent error measure, and define the 
adiabatic error to be 

,5[x(s)] = ||l/(x,)-V;d(xs)||. (10) 

Since \\{V — Kid)|V')ll !i 11^ ^ Kid||, where the norm on the 
right-hand side is the standard sup-operator norm (often de- 
noted II • lloo) im 



\\X\\= sup i/(i;|XtX|w) = maxcr,(X), (11) 

\v): \\\v)\\=l ^ 

where {(Tj(X)} are the singular values of X (eigenvalues of 
an upper bound on (5[x(s)] is then also an upper 
bound on \\\ij{xs)) - |?/'ad(xs))||. 
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Using the fact that the sup-operator norm is unitarily in- 
variant (||1/AVF|| = ||A|j for any operator A and any pair of 
unitaries V and W) we can rewrite S as 



5[x(s)] = ||/-n(x,)||, 

where the "wave operator" 
satisfies the Volterra equation 

n{s)^i- f KTis')n{s')ds', 

Jo 

with the kernel 



(12) 



(13) 



(14) 



Kt{s) = V^,{s)[dsPo{s),Po{s)]V,i{s). (15) 

Considering Eq. (j9|, —iKT{s)/T is simply the interaction- 
picture Hamiltonian which results from transforming H{s) 
to the interaction picture with respect to iJad(s), where 
i[dsPa{s), Pq{s)]/T plays the role of the perturbation. There- 
fore, in analogy to the Dyson series of time-dependent pertur- 
bation theory, the Volterra equation can be solved by iteration, 
which yields 



(16) 



1=0 



where 



Qois) = /, (17) 
^i>o{s) = - I KT{s')ni^i{s')ds'. (18) 



As shown in Refs. GUTtI, V/ G {2k - 1, 2k} [k € M): 

snv\\ni{s)\\^0{l/T^), (19) 

s 

l-l 

sup \\n{s) - J2 = 0{1/T''). (20) 

Using the above results, ||/— ri(s)|| can be expressed in terms 
of a 1/T series expansion, since 

\\i-n{s)\\ = \\n,{s) -Y, r KT{s')ni_,{s')ds'\\ 

r \\KT{s')\\Y,\\ni^,{s')\\ds' (21) 

•^0 i>2 



where 



|r!i(s)||+?(s)o(i/r), 



e{s)^ / \\[d,,Po{s'),Pois')]\\ds\ 
Jo 



(22) 



(23) 



Thus the error S is upper-bounded as 

<5[x(s)] <<5i(,s) + <52[x(s)], 

where 

6,{s) ^ \\n,{s)\\^0{l/T) 
S,[^{s)] ^ ?Ix(s)]0(l/T). 



(24) 



(25) 
(26) 



Both error components can evidently be made small by choos- 
ing a large T, while for a given T, 62 can additionally be made 
small by choosing a path over the control manifold A4 with 
small e. Note that in addition to ||rii(s)|| 0{l/T) we also 
have the bound < ||XT(s')l|ds' = elx(s)], but 

from Eqs. ^T9\ and pO] ) as such we do not have a bound of the 
form \\ni{s)\\ < ?{x(s)]C'(l/r). One can see fromRef. GSl 
how Si{s) depends on T, the gap, and the norm of the Hamil- 
tonian or its derivatives. However, the coefficient of the 1/T 
term of 5i does not appear to have a geometric significance, 
and we shall therefore exclude from our study of adiabatic 
geometry in this paper. 

In the following we shall make the upper bound on S small 
by finding a path which makes e|x(s)] small. Finding the path 
which minimizes e is, however, beyond the scope of this work. 
Instead, as we show below, after replacing the sup-operator 
norm by the Frobenius norm, the problem of minimizing 62 
has a geometric solution in the sense that a Riemannian metric 
tensor is encapsulated in e[x(s)] [Eq. (23 1 with the modified 
norm] . 

In Appendix |B] we prove that 



\[dsPo,Po] 



\Po{dsH) 



1 



H-Er 



idsH)Po\ 



(27) 

where [H — Eq]^^ is called the reduced resolvent and is a 
shorthand for (/ - Pq)[H - Eo]~\l - Po)- 

For a different method of traversing eigenstate paths of 
Hamiltonians, based on the use of evolution randomization 
and a quantum phase estimation algorithm, see Ref. I.20J . 



2. Nondegenerate case 

When H has a discrete and nondegenerate spectrum, Pq = 
|$o>(*o| and / - Po = E„>o l*n)($„|, where {|$„)}„>o 
are the excited eigenstates of H with eigenvalues {Pn}n>o- 
In this case 



1 



E 



1 



H — Er\ ^ — ^ En ^ En 

^ n>0 ^ ^ 



l*n)(1>n| 



(28) 



Using the chain rule of differentiation to write dgH = 
{diH)x\ where dot denotes dg and di denotes 9/9x*, and 
using the Einstein summation convention Eq. ( [27| is easily 
simplified in the nondegenerate case to yield: 

?lx(s)] - 72glj'(x)x'xJds', (29) 
Jo 
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where 



where the metric tensor is defined as 



^ {'i>o\diH\<^n){^n\djH\<^>o) 

.n>0 



(30) 



The manner in which g'j' appears in Eq. (29 1 suggests that it 
plays the role of a metric tensor. This metric tensor is identical 
to the metric tensor which was identified in the differential- 
geometric theory of QPTs 1 10|. We next consider how to gen- 
eralize this result to the degenerate case. 



Metric tensor for the degenerate case - moving to the 
Hilbert-Schmidt norm 



We would like to identify Eq. ( |27| i with a metric tensor 
However, the appearance of the sup-operator norm presents a 
problem, since this norm need not be differentiable. Hence 
we replace the operator norm with the Frobenius (or Hilbert- 
Schmidt) norm 



rank(^) 



which satisfies ifTSll 

\\X\\ < \\Xh < Vrank(X)||X|l. 



(31) 



(32) 



Note that the operator Pa{dsH) 



H^Eo 



{dsH)Pa appear- 
ing in Eq. p7| ) has support purely over the ground-state eigen- 
subspace of H, due to the projections Pq to the left and right. 
Therefore its rank is at most qq, and as a consequence of 



Eq. ( 32 1 the replacement of the operator norm by the Frobe- 
nius norm does not alter || [dgPo, Po]\\ (hence ?or g) for the 
nondegenerate case (go — 1), while it enables a differential 
geometric bound in the degenerate case, which is at most y/go 
times greater than the expression obtained with the operator 
norm. Additionally, and this is our main reason for moving to 
the Frobenius norm, it guarantees analyticity of the adiabatic 
error and the metric tensor when H is analytic. 

With these considerations in mind, let us now redefine the 
adiabatic error using the Frobenius norm: 

e{s)= f \\[ds'Po{s'),Po{s')]hds'. (33) 
Jo 

Then e{s) < e(s) < ^/go£{s) and consequently 

S2is) < e{s)0{l/T) < ^52{s). (34) 

Minimization of e(s) thus "squeezes" the error component 62- 
We show in Appendix [C| that 



^25og*j(x)x»xJds', (35) 



1 
2^ 
1 

2go 



Tr[d,PodjPQ] 



It is simple to verify that reduces to g^j'' in the nondegen- 
erate case, and similarly e(s) reduces to ?(s) in this case. 

Standard calculus of variations then tells us that minimiza- 
tion of e[x(s)] is tantamount to finding the geodesic path 
which is the solution to the following Euler-Lagrange (EL) 
equations: 



(36) 

+ i ^ j- 
(37) 



X* + r^fcX^x'^ = 0, 



where the connection T is 



(38) 



(39) 



We have thus endowed the control manifold M with a Rie- 
mannian structure, given by the metric tensor g : Ta^ (81 Tm 
R. That g really satisfies all the properties required of a met- 
ric is shown in Appendix [P] Other geometric functions such 
as the curvature tensor R can be calculated from g 11211 . 



C. Operator fidelity 

Another approach to the adiabatic error is provided by the 
"operator fidelity" ll22l between V and V^d, 



(40) 



where g is an arbitrary density matrix of the system, which 
here we take to be the totally mixed state I/N. The operator 
fidelity derives its name from the fact that it quantifies the fi- 
delity in the entire Hilbert space, and unlike our previous error 
measures e'and e, which involve the ground state projector Pq, 
is not restricted just to ground states. However neither is the 



adiabatic error 6 [Eq. ( 10 1] restricted just to ground states, and 



the two are obviously closely related. In Appendix|E]we show 
that 



1 



(41) 



so that minimizing e maximizes fg, and vice versa. 

Let O be an arbitrary observable, and consider it in the ro- 
tated bases associated with the actual or adiabatic dynamics: 

0(s) = V{s)OV^{s) (42) 
Oad(s) = Kd(s)0Kd(5)- (43) 
In addition to the bound (j4T]l we show in Appendix |E] that 

\\0(s) - 0.,,{s)\\ < \\0\\{S,{s) + 52[^{s)]) 

[2 + 0{l/T)], (44) 
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which is identical to the adiabatic error bound (24 1, apart from 
the factor \\0\\[2 + 0{l/T)]. Thus our bound of the operator 
distance || 0(s) — Oad(s)|| also has the component 5i and the 
component 82 with its apparent geometric contribution, which 
can be squeezed by choosing a geodesic path, as in subsection 
HlBll 



D. Natural geometric formulation 

1. Grassmannian 

An alternative, natural way to obtain a geometry for adi- 
abatic evolutions employs the Grassmannian structure of the 
dynamics |23|. As explained above, in the ideally adiabatic 
case the eigensubspaces corresponding to the ground state and 
the rest of the spectrum (Pq and / — Po> respectively) do not 
mix; each follows its own unitary dynamics determined by 
its Wilczek-Zee holonomy, hence V-^ = ® This 
implies a Grassmannian manifold 

GN.g, - U{N)/U{go) X U{N - go) 

= {Po e V{n) I Po' = ^0, Tr[Po] = go}, (45) 

where U{k) is the group of k x k unitary matrices, and ViT-L) 
is the convex space of all density operators (positive semidef- 
inite, unit trace matrices) defined over T-L. A natural distance 
(metric) over this space is given by 



d{Po,Po)^ 



0II27 



(46) 



whence, keeping only the lowest non-vanishing order, we 
have 



d2(Po(x),Po(x + dx)) 



1 

1 



|Po(x + dx)-Po(x) 



1 



— Tr[(dPo(x) + ^d2Po(x))2] 
1 

1 



Tr[dPo(x)dPo(x)] 
Tr[9,Podx'9jPodx^] 



2ffo 
gijdxMx-', 



(47) 



with the metric tensor as defined in Eq. (36 1. Thus the adia- 
batic metric tensor is precisely the metric over the Grassman- 
nian manifold defined by the ground state projectors. 



2. Adiabatic parallel transport 

In this subsection we wish to define a notion of adiabatic 
parallel transport. We start with the standard purification 11261 



|27l 



of Po, where U is an arbitrary unitary acting on T-L, so that 
Po = WW^ . Here W is considered a vector in a larger (ex- 
tended) Hilbert space Hcxt^ i e-. a pure state whose reduc- 
tion yields (the density matrix) Po. The Hilbert space Hoxt 
is equipped with the the Hilbert-Schmidt inner product 



{A,B) :=Tr[A^P]. 



(49) 



Given Po, the fiber of all purifications sitting on the unit 
sphere 5(Hext) {W e Hext : {W,W) = 1} of Hext is 
the Stiefel manifold of orthonormal go-frames of Hcxt^ where 
Tr[Po] = go (i.e., the set of ordered go-tuples of orthonormal 
vectors in Hcxt)- The gauge transformation (48 1 means that 



the fiber admits the unitaries of Ti. as right multipliers. Infor- 
mally, the Us act as arbitrary "phases" associated with Po. 

Starting with a curve of (unnormalized) density operators 
s !—> Po(s) and one of its purifications 



s^W{s) Pois) ^Wis)W\s), 



(50) 



the length iu[s] = J^{W{s'), W{s'))ds' of the curve in "Hcxt 
is not invariant against gauge transformations (48 1. The Euler 
equations for the variational problem £[s] mfjj £u[s], i.e., 
for the geodesic are 1 26 , 27 1 



(51) 



also known as the Uhlmann parallel transport condition. Sub- 
stituting W^t ^ [/tP(j and dW = (dPo)C/ -I- PodU yields the 
condition 

U^Po{idPo)U + PodU) = {U^dPo + {dU^)Po)PoU, (52) 

which, using UdW = —{dU)W, reduces to 

PoidU)U^ + idU)U^Po - [dPo, Po] (53) 

on the vector bundle over the Grassmannian Q jv.gn . Here U = 
U{s) is a general unitary undergoing parallel transport as s i-> 
Po(s). We now seek those unitaries U which in addition to 
parallel transport, also satisfy adiabaticity. 

To this end let J(s) be the infinitesimal generator of U (s), 
i.e., 

idsU{s) =TJ{s)U{s). (54) 
Substituting this expression into Eq. ( [53| we obtain 

(55) 
(56) 



PoJ + JPo = i[dsPo,Po]/T 

= Pt-dd — H, 



where in the second line we used Eq. (|9|. Thus, U satisfies 
adiabatic parallel transport if in addition to being a solution 
to the parallel transport condition ( [53] l its generator also satis- 
fies the adiabaticity condition 

Po J + JPo = 0. (57) 

What is the generator J which satisfies this last condition? 



W = PoU 



(48) Using Eqs. (BlOi and (B12i for the nondegenerate case we 
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obtain 



tT{PoJ + JPo) = 



[Po,Po] 

1 



E 

n>0 



H-Eo 

Pom 



HP„ + P„H 



1 



H-Eo 



En — E( 







(58) 



Taking matrix elements we find ($ol>^l*i'o) = and 
-ir($o|J|$fc> = ;g-^($o|i^|*fe), while the matrix ele- 



ments of J between the excited states are unspecified, so that 



T ^ En- Eo 



|$0)(*n|+H.C. + Ji, (59) 



where J± is an arbitrary operator satisfying J± = QoJ±Qa- 
Instead of trying to obtain perfect adiabaticity (_ffad = H) 



we can settle for an approximation. Noting that Eqs. (33 1 and 



(55 I imply 



e{s) = T / \\PoisW) + J{s')Pois')\\2ds', 



(60) 



it follows that minimizing e, or equivalently finding the adia- 
batic geodesic, endows the "phase" U of Pq with an adiabatic 
characteristic which is compatible with the Uhlmann paral- 
lel transport condition. Thus, we have shown that the metric 
tensor g emerges naturally also from the notion of adiabatic 
parallel transport. 



3. Bures metric 

There is also a straightforward connection between our 
metric and the Bures metric f28l|. For two arbitrary density 
matrices pi and p2, the Bures distance is defined as 



rf2 



s(Pl,P2) = 1 - F[pi,p2), 



(61) 



where F{pi,p2) = ^t\{p\^'^ p2p]!'^fl'^\ is the fidelity be- 
tween these two states ll29l [30l . When the density ma- 
trices depend on a parameter x, the infinitesimal distance 
'^Bures (p(x), p(x + dx)) Can be shown to be ||4J 



,(p(x),p(x + dx))-Tr[p(x)£2(x)], 



(62) 



where £(x) is the "symmetric logarithmic derivative," (SLD) 
defined via 

dp(x) = i (/:(x)p(x) + p(x)£(x)) . (63) 

From the property Pq = Pq we obtain 

dPo(x) = dPo(x)Po(x) + Po(x)dPo(x). (64) 



and hence [see Eq. ( A2 1] 

dgo = Tr[dPo] = 2Tr[PodPo] = 2Tr[PodPoPo] = 0, (65) 

i.e., the degeneracy is constant. Thus if p(x) = Po(x)/go, 
then d[Po(x)/5o] = Po(x)/.godPo + dPoPo(x)/5o, and the 
definition of the SLD [Eq. (|55]l] yields 

/:(x) = 2dPo(x). (66) 

Inserting this back into Eq. (|62]) results in 



4u>-es(^o(x),Po(x + dx)) = -Tr[Po(x)(dPo(x))^] 

= gf/-^(x)dxMx^ (67) 



where 



g^^"-(x) = -Tr[a.Po(x)9,Po(x)]. 
.90 

By comparison with Eq. (|36]), we obtain 



Bures 

5ij 



(68) 



(69) 



We note that the Bures metric is also connected to "quan- 
tum Fisher information tensor," which plays a principal role 
in quantum estimation theory [4]|28][3T1[32 1. In fact the Bures 
metric is (up to an unimportant constant multiplicative factor) 
equal to the Fisher tensor. Therefore, the adiabatic metric is 
the quantum Fisher metric. The role of the metric g in quan- 
tum estimation theory is thus highlighted naturally this way. 



E. Comparison of adiabatic metrics 

In adiabatic evolution (as well as in adiabatic quantum com- 
putation) 8 and T are the primary objects of interest. Our 
method for obtaining the metric g here is based on minimizing 
an upper bound on the adiabatic error 8 for a given evolution 
time T. In Ref. [11] we pursued a complementary route and 
proposed a different metric. 



g.,(x) = Tr[a,i/(x)a,i/(x)]/A4(x) 



(70) 



derived from minimizing a time functional inspired by the tra- 
ditional adiabatic condition. We called this the "quantum adi- 
abatic brachistochrone". 

The major difference between these two metrics is in their 
distinct gap dependence. This can be understood, for example, 
by noting that 



whereas 



< 



< 



min, 



\d^Hd,H\\r 
min, 



(71) 



(72) 



where ||X||i = Tr[\/XtX] = Yl,i^iW is the trace norm 
|[T8l| (see Appendix [G| for the proof). Thus, the metric g has 
a quadratically less dependence on the inverse gap. This may 
imply different behaviors for these metrics and their corre- 
sponding curvatures; hence they are essentially distinct. 
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F. Strategies for reducing tlie adiabatic error and their effect 
on geometry 

Considering that g is related to minimizing the upper bound 
on 6, it is useful to briefly recall how 6 scales with T and how 
this scaling may be improved. 

Rigorous proofs of the adiabatic theorem — based on suc- 
cessive integration by parts of O — state that if {H{s)} is 
a family of {k times continuously differentiable) in- 
terpolations/paths with bounded {I E {l,...,fc}) 
and compactly supported dgH over s E (0, 1), then 6 — 
^(l/r^C^-i)) |[l6l[l7l[l3- If these assumptions are supple- 
mented with that of analyticity of H{s e C) in a small strip 
around the real time axis, and if in addition 



diH{0) = diH{l) =Oyi<k, 
the result is an exponentially smaller error 

6 = 0{e-^^), 



(73) 



(74) 



where c = min^ A^/ max^ ||9s-H^|P (up to an 0(1) prefactor) 
I33J. 

Our path — as the solution to the second-order differen- 
tial equation ( [38] l — minimizes e rather than 6, which is 
not necessarily compatible with the boundary conditions 
dgH{{Q, 1}) = 0. Thus, in principle, there remains room for 
further optimization of the path for S beyond what is captured 
by simply minimizing its upper bound e{s)0{l/T) |16, 17|. 
Such finer optimizations, however, may not always result in a 
Riemannian geometry because the corresponding functionals 
and Euler-Lagrange equations would depend on higher deriva- 
tives of H. 



QPTs take place at zero temperature, in which the system is in 
its ground state, quantum criticality should be identifiable by 
ground state properties. Specifically, the ground states right 
before and right after a quantum critical point are expected to 
have very little overlap. In this manner ground-state fidelity 
may be considered as a natural, fairly general order parameter 
for quantum critical systems, irrespectively of their internal 
symmetries lilOJ . We shall discuss this feature in more detail 
below. 



A. Metric tensor for QPTs 

Here we derive the metric attributed to QPTs for the case of 
degenerate ground states as a natural extension of the similar 
metric proposed for the nondegenerate case 1 10 1. 

In the degenerate case we should work with the ground- 
state projector Pq{x.). A variation in the properties of Po{x), 
caused by the change x — > x + dx in the Hamiltonian pa- 
rameters, can be captured by the order parameter chosen to be 
the operator fidelity of Pq (x) and Pq (x + dx) relative to, e.g., 
I go (J go is the go X go identity matrix), 

/,(Po(x),Po(x + dx)) - (Po(x),Po(x + dx)), 

= 1 - G,j (x)dxMx^ (75) 



in which the Hermitian matrix 



= -Tr[Po(9^Po)(a,Po)^o] 
5o 



(76) 



is the "geometric tensor" for the degenerate case (see Ap- 
pendix [F| for the proof). Thus the information about the crit- 
icality of the quantum system is contained in the G tensor 
III. CONNECTION TO QUANTUM PHASE TRANSITIONS Note that in the nondegenerate case (go = 1) G.y reduces to 



The other physically-important aspect of our geometric for- 
mulation emerges from the observation that the metric g also 
arises naturally as the underlying geometry of QPTs. QPTs 
take place at zero temperature |7 j, where the system is in prin- 
ciple in its ground state. Such phase transitions exhibit pecu- 
liar behaviors and "orders," radically different from their ther- 
mal counterparts. In particular, in contrast to thermal phase 
transitions, the standard paradigm of the Landau-Ginzburg 
symmetry -breaking mechanism fTl l??! fails to explain the un- 
derlying physics of some QPTs. In fact, defining an appro- 
priate local "order parameter" — an essential ingredient of the 
Landau-Ginzburg theory — is not straightforward for a quan- 
tum critical system; some QPTs, such as those involving 
"topological order," provably do not admit any local order pa- 
rameter 1351 [36l . Additionally, tracking singularities of the 
ground-state energy cannot always foreshadow QPTs; quan- 
tum systems with matrix-product states may elude this test 

Gil. 

Notwithstanding the above subtleties with identifying 
QPTs, it has recently been shown that the simple notion of 
the "ground-state fidelity" is remarkably successful in signal- 
ing QPTs Iil0ii38i . This can be understood by noting that since 



G„- = (5i$o|5j$,)-(a,$,|$o>(*o|9j$o). (77) 

Accordingly, a Riemannian QPT metric tensor can be de- 
fined through 



gQ"(x) ^ Re[G.,(x)] 



1 

Wo 



Tr[5,PoajPo] 



I], 



(78) 



wher e w e used the same trick as that used in arriving at 
Eq. (C2i. Therefore, the QPT metric tensor g'^''^ is the 



same as the adiabatic quantum evolution metric g defined in 



Eq. (36 1 



B. Quantum critical scaling of the QPT metric tensor 

The critical behavior of a quantum system with a degen- 
erate ground state can be characterized by the metric tensor 
g. This is already evident from the fact that the divergence 
of g'5^^ is a sufficient condition for signaling a quantum crit- 
icality. To further elaborate on this connection, we follow 
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Ref. 113911 and obtain the scaling of the geometric tensor (|76 

2 



1 

go 



-Tr 



1 



H - E, 



-, 90,97, 

90 ^ ^ 



n>0 a,?7— 1 



(79) 



,(80) 



and via Eq. (78 i also for g^ . For simplicity we restrict our- 



selves only to gapped quantum systems with second-order 
QPTs. Thus, in a critical region x « x^,, the correlation length 
^ and the gap A exhibit the following scalings 



(81) 



with the critical exponents v > and zu, where z > is the 
dynamical exponent Q. The geometric tensor G has an inte- 
gral representation which not only facilitates the derivation of 
the scaling relation for G, but also enables an interpretation 
for G in terms of correlation (or response) functions. Indeed, 
as shown in AppendixlH] Eq. dSOll can be expressed as 



(82) 



Gij = — / dTTe-'P-'(TT[Pod^HrdjH] 
90 Jo ^ 



Tr[Pod,H]Tv[Pod,H] 

90 



p=0 



with d^Hr = e'^^d^He-^". 

Now we make some generic assumptions about the Hamil- 
tonian H. First, let diH he a local operator; that is, one can 
write 



(83) 



in which y labels the spatial region over which the local op- 
erator hi{y) has support. Second, the hi{y) operators have 
well-defined scaling dimensions ai near the quantum critical 
point Xc, such that if 



y ^ ay, r o^r, 

for a > 0, we obtain 

h^{y) a-"^h,{y). 



(84) 



(85) 



Under these transformations, Eq. ( [82[ ) yields the following 
scaling for the rescaled geometric tensor in the thermody- 
namic limit 



where 



Kij — ai + aj — 2z — d. 



(86) 



(87) 



Here, L is the linear size of the system and d is its spatial 
dimension. From Eq. (81 1, we obtain |x — x^ ^ C^^^"; 



i.e., the scaling dimension of the Hamiltonian parameter x is 



l/iy. Following standard scaling analysis arguments, the scal- 
ing behavior of the metric tensor (recall that g — Re[G]) in 
the off-critical limit ^ L is 



L-^lx- 



(88) 



Moreover, in the critical region, where £^ ^ L ^ the spacing 
between adjacent particles on the system lattice, in addition to 
the regular extensive scaling L"^, the finite-size scaling of the 
metric is g^ ^ ^rf-tij^ which could be extensive, subexten- 
sive, or superextensive [k ~ 0, positive, or negative, respec- 
tively]. We also remark that there exist models, exhibiting 
quantum topological order, in which the critical g scales log- 
arithmically, e.g., g ^ In |x — Xc| li4&H42l . 



IV. ADIABATIC GEODESICS 



In this section we solve the geodesic equation (38 1 an- 
alytically for some specific examples. Note that since the 
eigenprojections do not depend on Tr[i/], Eq. (38 i corre- 



sponds to an underdetermined system of coupled second- 
order differential equations. This can be seen more clearly 
by adopting a new parametrization (i.e., coordinate system) 
y(x) for the Hamiltonian such that ff (y(x)) = y^(x)]l + 
ff'(y2(x), . . . ,y^^(x)), in which y^ = Tt[H]/N and H' = 
H - Ti-[H]t/N. Since Fo(y) does not depend on y\ the 
metric g(y) does not depend on this parameter either. Inde- 
pendence from y^ translates in terms of x into the statement 
that only M — 1 equations in the system ([38]l are independent. 



A. Examples 

1. Deiitsch-Jozsa algorithm 



In the Deutsch-Jozsa algorithm B3l one is given an oracle 
that calculates a function / : {0, 1}" i— >■ {0, 1}. The promise 
is that / is either "constant" or "balanced," meaning respec- 
tively that, /(!) = /(i') Vi,i' or /(half of alli's) = [291 
The objective is to conclude whether / is constant or balanced. 
The Deutsch-Jozsa algorithm finds the answer by querying the 
oracle only once, while classical deterministic algorithms re- 
quire a number of queries that is exponential in n. 

An adiabatic version of this algorithm was introduced in 
Ref I 44J . We consider the unitary interpolation Hamiltonian 

m 



if(x(s)) =F(x(s))Fo0(x(5)), 



(89) 



where t/(x(s)) = 6*2'^'^''^'^, in which the Hermitian/unitary G 
operator is defined by G|i) = (— 1)-'^^'^ Here Hq is chosen 
such that |$o(0)) = 1+)*^" = 2-"/^ J2f2~^ |i) is its ground 
state, e.g.. 



^0 - ho J2k= 



(90) 



where |±) = (|0) ± |1))/V2, a, = |0)(0| - is a 

Pauli matrix, and /iq > is an energy scale. The bound- 
ary conditions are chosen as (xo,xi) = (0,1) such that 
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H{0) = Hq and H{1) = GHqG^; the latter guarantees that 
|$o(l)) = G|$o(0)) is the ground state of H{1)^ 

From Eq. ([89l it is seen that |$o(s)) = V^(s)|$o(0)), 
whence we obtain 

Po(x(s)) = 2-«EM'=oe^*'^(^)[(-^)'*"-(-^^'"''l|i)(i'|, 

axPo(x(s)) =^f[G,Po(x(s))]. (91) 

A straightforward calculation then yields 

g = Tr[(a,Po(x))'] 

= ^(Tr[Po(x)GVTr[(Fo(x)G)^ 

2"-l 

l-2-2"(^ e"^W)" . (92) 



vr 



i=0 



Since g is independent of x(s), the geodesic equation (38 i 
reduces to x = 0, whence the geodesic is simply 



x(.) 



(93) 



which corresponds to a rotation of the initial Hamiltonian Hq 
at a constant rate. 



2. Projective Hainiltonians 



Consider the following Hamiltonian: 

77(x(s)) = xi(s)P,^ +x2(s)P^, (94) 

where P^^ = 1 — |a)(a|, for a given |a) E H (similarly 
for Pj^), (a|b) is a given function of N, and the boundary 
conditions are xq = (1,0) and Xi (0, 1). This Hamilto- 
nian may represent the adiabatic preparation of an unknown 
("hard") state |b) from the supposedly known ("simple") ini- 
tialization I a) , provided that one has access to the "oracle" 
ifTTl . An important instance of this class is Grover's Hamilto- 
nian for search of a "marked" item among N unsorted items 
PI-6J (generalized to arbitrary initial amplitude distributions in 
Refs. illllll), where |a) = J2k=o \k)/y^ and |b) = |m), 
for TO G {0, . . . , TV — 1}. A successful adiabatic version of 
this algorithm was first described in Ref . [49J . 



Expanding Eq. (94 1 in the {|a), |a^)}^j^^ basis and using 
Eq. 05\ yields 



ff(x) 



x2(l - laoH -x^aoai 
— x^aocki x-'-+x^|ao| 

(x^ +X^)/{2,...,Ar_i}, 



(96) 



where we have used the completeness of the basis to write / — 

|a)(a| + Ia^)(a^|, I{2,...,n-i} ^ ^.=2' |a^)(a^|, 

and the matrix on the right hand side is written in the 
{|a), |a;j'-)} (sub-)basis. It then follows from Eq. (96i that the 
spectrum of H consists of the two nondegenerate eigenvalues 



E± 



1 



(xi+x2 



± V(xi)' + (x2)2 + 2(2|(a|b)|2 - l)xix2), (97) 
and the {N — 2)-fold degenerate eigenvalue 



^x^ 



(98) 



Thus, the gap between the ground state (P-) and the first ex- 
cited state (P+) becomes 

A(x) = ^{x^y + (x2)2 + 2(2|(a|b)|2 - l)xix2. (99) 



The Hamiltonian ( 96 1 can be diagonalized by noting that 
one can rewrite 



1 



i?(x) = -A(x)[A(x)I],-(xi+x2)/{o,ij]At(x) 



(100) 



+ (xl+x2)/, 



where is the Pauli matrix — diag(l, —1) e 
|1)(1| padded with zeros to embed it trivially into the A^- 
dimensional representation (i.e., = diag(o'2, 0, . . . , 0)), 
I{o,i} = diag(l, 1, , 0, . . . , 0), and the 2x2 unitary matrix 
yl(x) is defined as 



A(x) 



-i9{yi)(7y 



(101) 



(the extension to N dimensions is similar to that of by 
padding with sufficiently many zeros) with 



Since the Hamiltonian (|94]l is effectively two-dimensional 0036' = 2x^|(a|b)|v/l - |(a|b)|2/{4|(a|b)|2(l - |(a|b)|2) 



over the span of the vectors |a) and |b), it can be diagonalized 
analytically. Indeed, given |a), we have the freedom to choose 
A^ — 1 vectors {|a^)}^^^ at will such that together with |a) 
they constitute an orthonormal basis for T-L. I.e., (a|a^) — 
and (a^|aj-) = 6ij. Thus we can decompose |b) = ao|a) + 

^fS^ Qfija^). Utilizing the freedom in choosing {|a^)} (up 
to the orthonormality condition), we can always rotate them 
such that aiyi — 0. In this case, we have 



(x^)2 + (xi-[l-2|(a|b)|V-A)^} 



2-, -1/2 



(102) 



After removing the energy shift (x^ + x2)/ from Eq. ( 100 1, it 
is evident that the ground-state projection is 



(95) 



where ao = (a|b) and ai — (a]'-|b) (or more explicitly: 
ai = e*"^! ^1 - |(a|b)|2, for some arbitrary e [0, 2tt)). 



Po(x)=A(x)|l)(l|At(^)^ 
(padded with zeros). This yields 

g,, = Tr[a,PoajPo] 

= d^edjO. 



(103) 



(104) 
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Obtaining the geodesic for the one-dimensional case x = 
(1 — X, x) turns out to be simple and can be performed analyt- 
ically, yielding 

x(s) = - , l^'^*^^* tan[(l - 2s) arccos I (alb) 11. 

^ ' 2 2^/1- |(a|b)|2 ' ' 

(105) 

It is interesting to note that this is exactly the solution obtained 
in Ref. ifTTl from the different metric g [Eq. (70 1]. 



To make further progress we focus on the one-parameter 
cases: (i) x = (1 — x, x), (ii) x = (x, 1), and (iii) x ~ (1, x), 
all subject to the boundary conditions x(0) = 1 — x(l) = 0. 

Let 



p(x) 



1 



4^[1-2(1 



2lTi 



2m+ 



T)(i-x)xr 



(111) 



For a given finite lattice size m, the geodesic equation for case 
(i) reads 



3. One-dimensional transverse-field Ising chain 

Consider a one-dimensional chain of spin-1/2 particles in- 
teracting according to the following Hamiltonian: 

m 

i/(x(s))^- Y: xi(s)aW+x^(s)aWai^+i), (106) 

£— — rn 

with the boundary conditions xq = (1, 0), xi — (0, 1), and 
^(m+i) ^ ^(1) |fjQ| £xact diagonalization by the Jordan- 
Wigner transformation |7| yields 



2p(x)X + 9,p(x)(x)2 = 0. 
This equation can be integrated to yield 



(112) 



2s 



x(s) .1/2 

v/R^dx'/ / v^K^dx', (113) 
Jq 



We next consider the thermodynamic limit m — > oo, where 
we can obtain a simple closed-form formula for the geodesic. 
The expression in this limit follows from substituting J^e ^ 
dz [with Zi = 27r£/(2TO + 1)] and taking into ac- 
count that the model exhibits a QPT at Xc = 1/2 correspond- 
ing to Sc = 1/2. This yields 



I $0 (x)) = ®T=i (cos Oe (x) 1 0) |0), + z sin 0, (x) 1 1) 1 1),) , 



(107) 



i(l-tan2[f(l-2s)]), 0<s< i, 
i(l+tan2[|(l-2s)]), i < s< 1. 



(114) 



where (cf Ref. flOl) 



sin 29£ 



x2 sin(^) 

^ 2m+l / 



J(x2 COS ^ - Xl)2 + (x2)2 sin2 ^ 

Y ^ 2m+l / ' \ / 2m+l 



.(108) 



For details of the derivation see Appendix |IX| 

Similarly, for both cases (ii) and (iii) we obtain the geodesic 
for a given finite m as 



It is evident from Eq. ( 107 1 that 

2 



K{S) 



v/^dx7 / v^dx', 



where 



1=1 

2 



4 ^ [1 

i=i L 



sin2( 2M^) 

~2i^t 12 • 



2 COS 



(115) 



(116) 



2m+lJ 



= 5Z 51 ^^^^ 1 0) + ^ COS 0, 1 1) 1 1) ,) 



1=1 t=\ 



In the thermodynamic limit a quantum critical point emerges 
at Xc = 1 (sc = 1), and a similar approach as in case (i) yields 
the geodesic 



where is the same as |<I>o) [Eq. (107i] except that the term 



x(s) = sin(7rs/2). 



(117) 



with the label I is absent. In addition it is easily verified that 
(<i>o|$o) = 0. Thus, we obtain 

2 m 



For details of the derivation again see Appendix IX 



After inserting these results into Eq. ( 36 1 we have 

m 

g,,(x) =^9,0,(x)a,0,(x). (110) 



This, then, is the geometric tensor for the transverse field Ising 
model. 



Figure[T]illustrates the geodesies obtained for the transverse 
field Ising model subject to the three parametrizations we have 
discussed. 



B. Geodesic for passage through a quantum critical point 

A limitation of our formalism is that, in principle, exact 
knowledge of the ground state is required in order to ob- 
tain the geodesic. Unfortunately, such knowledge is rarely 
available, the exceptions being certain exactly solvable mod- 
els such as those we treated in the previous subsection. With 
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0.0 0.2 



0.4 0.6 0.8 
s 




FIG. 1; (Color online) Optimal adiabatic paths for the one- 
dimensional transverse-field Ising model, corresponding to the pa- 
rameterizations x = (1 — x, x) (left), x = (x, 1) and x = (l,x) 
(right). The red dashed lines represent the thermodynamic limit, 
while the solid blue lines correspond to m = 1, 4, 10, 30, 100, ap- 
proaching the dashed line as m increases. 



partial knowledge or an approximation for the gap, one should 



solve Eq. ( 38 i on a case by case basis, possibly numerically. 

However, while these observations apply in a setting where 
one wishes to obtain the geodesic over the entire parame- 
ter manifold, the situation in the vicinity of a quantum crit- 
ical point is rather different. Indeed, the most interesting 
physics usually happens in the vicinity of the quantum crit- 
icality. In addition, the behavior of a quantum adiabatic algo- 
rithm is essentially governed by how the system approaches 
and/or passes through a quantum critical region. These con- 
siderations suggest that knowledge of the geodesic around the 
quantum critical region should suffice for most algorithmic or 
physically relevant applications, thus obviating the need for 
knowing Pq everywhere. 

Computation of the critical behavior of other geometric 
functions, such as F and R, is straightforward. E.g., in 
the one-parameter case, where x = (x), the Euler-Lagrange 



(geodesic) equation ( 38 i in the critical region slightly before 
and after the critical point reduces to x + i^kx^/2x — 0, 
whence 



x(s w Sc) ~ Xc + A{s - 



(118) 



After using a ~ d + z — l/v fST), where a is the scaling 
dimension [recall Eq. ([85|], we obtain 



(119) 



X = 2/(2 + i/k) = 2/diy > 0, 



with A constant (derivation details are given in Appendix|IX|i. 
This is a remarkable result as it characterizes the optimal adi- 
abatic passage through a quantum critical point in terms of 
the universality class of the system. Moreover, this result con- 
firms that the critical geodesic has a power-law dependence on 
s, as first reported in Ref. |52|, although away from the crit- 
ical region the dependence can be different. References ll52l - 
l55l report critical behaviors of the metric tensor and related 
parameters obtained using different methods, such as mini- 
mizing exact expressions for transition probability in thermo- 
dynamic limit. In contrast to the result of Ref. |52|, in our 
analysis the exponent x of the critical geodesic depends on the 
dimensionality d, whereas it is independent of the total time 
T. In adiabatic evolution the dependence on T is of course 
expected; however, note that our scaling result depends only 



upon the geometry of the control manifold, which does not 
depend on T. 



V. SUMMARY AND CONCLUSIONS 

In this work we set out to elucidate the role of geometry 
in adiabatic quantum evolution. By splitting the "adiabatic 
error", i.e., the norm of the difference between the ideal adi- 
abatic evolution operator and the actual propagator, into two 
components, one of which is endowed with a geometiic mean- 
ing, we were able to derive a Riemannian metric tensor which 
encodes the geometry of adiabatic evolution. This metric is 
capable of describing evolution over both nondegenerate and 
degenerate subspaces. We then showed that this same metric 
tensor arises naturally from a number of different but comple- 
mentary viewpoints, including a minimization of the operator 
fidelity, and a focus on the Grassmannian structure of the dy- 
namics. 

Our second major goal in this work was to establish a firm 
connection between adiabatic evolution and quantum phase 
transitions. By analyzing the infinitesimal variation in the op- 
erator fidelity we showed that, in fact the same metric tensor, 
arises in both cases. We further derived the quantum critical 
scaling of this metric tensor 

Having established a unified geometric framework for adi- 
abatic quantum evolution and quantum phase transitions, we 
proceeded to find the geodesies on the manifold described by 
the unifying Riemannian metric tensor Such geodesies are 
of particular interest in adiabatic quantum computing, where 
they correspond to paths which minimize the geometric com- 
ponent of the deviation between the actual and desired final 
states. We analytically determined the geodesies in three ex- 
amples of interest: the Deutsch-Jozsa algorithm, a general- 
ization of Grover's algorithm, and a model described by the 
transverse field Ising model. While such examples are im- 
portant as proofs of principle, one cannot in general hope to 
analytically find the geodesies. For this reason we focused on 
the passage through the quantum critical point, and showed 
that in general, for second order QPTs, the geodesic in this 
case obeys a universal scaling relation. 

Among other applications, we expect that the formalism we 
have developed will lead to further developments in adiabatic 
quantum computing, where the role of criticality is well appre- 
ciated. We expect additional applications in holonomic quan- 
tum computing, where degeneracy plays an essential role, and 
where a differential geometiic analysis of gate error minimiza- 
tion has not yet been carried out. 
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Note added. — While this work was being finalized for sub- proposes a generalized quantum geometric tensor related to 
mission, a related manuscript appeared |56 |, which similarly adiabatic evolution of quantum many-body systems. 



Appendix A: Proof of the Wilczek-Zee holonomy formula 

Notice that from the fact that Pq{s) is a projector, i.e., P(){s) — Pq{s), we obtain 

Pois) = Po(s)Po(s) + Po{s)Pois), (Al) 

(where Po{s) = dsPo{s)), so that 

Pois) Pois) Pois) = 0, (A2) 

and 

[Pois), Pois)] = 2Po(s)Po(s) - Ms). (A3) 
Let Qois) denote the projector orthogonal to Pois), i.e., Pois) + Qois) = /. Then we have 

Pais)Qois) = Qois)Pois) = 0. (A4) 
The differential equation for vj^^, (s) [Eq. ^] can be obtained as follows: 

dsVilis) = {^^is)\V.,is)\<i>^'iO)) + (<i>o"(s)|V;d(s)|<(0)). (A5) 
In addition, consider the action of i?ad(s) [Eq. (j9]l] on |$q (s)), 

HM\m^)) = {His) + 2iPois)Pois)/T-tPois)/T)\^^^is)) 

= Eois)\<^^is))+tPois)\^^is))/T. (A6) 
Since Pois) = E^li \H is)) {<i>^o i^)\ + |<f(^(s))(<i>^(s)|, we have 

90 

Pois)mis)) = \Kis)) + J2i^^o{s)\^oim<is)). (A7) 

0=1 

Using Eq. ( A6 1 and ( |A7[ ), we can rewrite Eq. ( |A5| l as 

dsvilis) = (<i>^(.)|v;d(.)|<(o))-zr($^(.s)|i/.d(s)v;d(.)|<(o)) 

go 

= -iTEois){-^>^is)\V,,is)\<^>^\0))~Y.(^^(')\^ois)){His^^^^^ (A8) 

/3=1 

Without loss of generality, after setting £'o(s) — 0, we obtain the following differential equation for V^^, (s): 



d.vilis) = -J2ims)\His)){<is)\K,is)\n'iO)) 

0=1 

90 

= (A9) 



13=1 

whose solution is 



with 

^/3 



(AlO) 



is) = V exp Ais')ds'^ , 

A^^ = i<i>^\ds\<^^o)- (All) 
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Appendix B: Proof of Eq. j27) 



Equation ( A2 1 yields 



[Po,Po] 



(PoPoPo + PoPoPo)- 



(Bl) 



Using Eq. ( Al i to write PqPq = Pq — PqPq and substituting this into the first term of Eq. (Bl 



we then have 



[Po,Po 



-{P^ -P^Po + PoP^Po) 
-Pq + QaPoPo- 



(B2) 



The second term vanishes, as can be seen by using Eq. ( Al i to write Pq = {PqPo + PqPo)^'- 

QoPqPo ~ Qo{PoPoPoPo + ^o-Pq^o + PoPqPq + PoPoPoPo)Po = 0, 



where we used Eq. ( A2 1 on the first two summands and Eq. ( A4 1 on the last two. Thus we conclude that 

[P,,Po]' = -Pl 



(B3) 



Note that Pq = X^aLi l*o)(*ol + l*o)(*ol is Hermitian and that therefore [P(),Po] is anti-Hermitian. Thus both Pq 
and [-Poi^o] are unitarily diagonaUzable: —Pq — VDV\ [Pq,Pq\ = W EW \ where V and W are unitary, while D and E 
are the diagonal matrices of eigenvalues. Therefore it follows from Eq. Q that ||FD2yt|| ^ || VK£;2W^^||, and from the 
unitary invariance of the operator norm that 
eigenvalues are equal, i.e.: 



From here we conclude that the maximum absolute values of their 



It also follows that HPq^H = || [Po, ^o]^ 



Next we wish to show that 



m.p,]\\^m\ 
\Dr = m,p^]r 

m,p^]\\ = Jiii\ 



, I.e., 



Pn = - PuH 



1 



1 



-HPr 



(B4) 



(B5) 



(B6) 



M — Pq P^ — Eq 
To prove this note first that the Hamiltonian can be decomposed as 

H = EaPo + QoHQo. (B7) 

Then 

H = EaPo + EoPo - PoHQo + QaHQo - QoHPo, (B8) 
and multiplying this equation by Pq from the right while using Eqs. (|A2|i and (|A4[) and the fact that H commutes with Pq, yields 



HP„ = EoPo + EoPoPo-{I-Po)HPnPo 
= EoPo + EoPoPo-HPoPo. 



(B9) 



The operator H — EqI is invertible when its domain excludes the spectrum of H (and is then called the "reduced resolvent;" 
see, e.g., Ref. Il33l ). That is, the inverse is defined as Qo[Pt ~ -E-ol^Qo (but for brevity and when there is no risk of confusion, 
we simply write [H — Eq]^^ henceforth). With this restriction in mind we then have 



PoPo = -- 



H- E, 



^ [H - Eo)Po = - ^ 



where in the last step we used 



H -E, 



H-En 



H - E, 



0, 



HPo, 



(BIO) 



(BID 
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which is due to the fact that the range of [H — Ea] ^ is the range of Qq [recall also Eq. (28 i]. Similarly, by multiplying 
Eq. (|B8[) from the left by Pq we obtain: 



PoPo = -PoH 



H-Eo 



(B12) 



Adding Eqs. (BlOi and (B12 1, and using Eq. ( |A1[ ) again then yields Eq. (B6i 
As a corollary, we can also calculate E(){s) from Eq. (B9 1 



Eo{.s)^TT[HPo]/go. (B13) 
Calculation of Pq or higher order derivatives of Pq follows similar logic (see, for example, Ref. 1,33] ). For example, we obtain 



1 

H-Er, 



1 

H-En 



HPo + ^ HPn + ^ HPn 

11 — JZ/Q 11 — £/Q 



This relation can be simplified further after replacing Pq [Eq. (B6 1], using the identity 



1 

H-Eo 



^- (H -Ea) ^- , 



(B14) 



(B15) 



and inserting Eq [Eq. (B13 i]. However, we do not need the final explicit form here. 

1 



We are now ready to prove Eq. (27 1. Let 



A = 



1 



H-Eq 



HPq, B = PqH 



H-Eq 



(B16) 



Then, using Eqs. (B5 i, (B6 1, and (Bll i yields 



\[Pa.Po 



(B17) 



Note that A and B are both positive operators and that they have orthogonal support. Therefore + = 

max{ 1 1 A 1 1 , 1 1 i?^^ S 1 1 } . Moreover, we have A^ A = BB\ and it is a ba sic p roperty of the operator norm that = 



for any operator B. Thus yJ\\A'iA + B^B\\ = yJ\\A^A\\, which is Eq. (27 1 



Appendix C: Proof of the error formula in tlie Frobenius norm 



Starting from the definition of the adiabatic error, Eq. (33 i, we have, by using Eq. (A2i together with Pq = Pq and cyclic 
invariance of the trace: 



e(s) = / \/Tr[{PQPQ-PQPQ){PQPQ-PQPQ)]ds' 

Ti[PqPqPq + PqPqPoW 

TT[Po{d,PQ)id,PQ) + id,PQ)PQ{djPQ)]±^±lds' 

Jo " 

where Pq ~ diPox.^. Using Pq — Pq once more to obtain -Po(9iPo) + {diPo)PQ ~ OiPq we have: 
e[s) 



y^Tr[{5,Po - {d,PQ)PQ}{d,PQ) + (9,Po)i'o(9jPo)]x»i^ds' 
\/2.gogy(x)x*xJds' 



(CI) 

(C2) 
(C3) 



where the metric tensor is defined as g^j = Tv[diPQdj Pq] / 2gQ, which is Eq. (36i 
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Next let us derive Eq. ( 37 1. From Eq. (|B6|l we have 



a,Po - - Poi^^H) 



H — E(:\ H — E( 



{^^H)Po 



Inserting this into TT[diPodjPo] and expanding the product while using Eq. (Bll 

TT[d,Pod,Po] = Tr 



Tr 



Tr 



Poi^^H) 



1 



1 



H — Eq H - Eq 
1 \ / 1 



we obtain: 

{d,H)PQ ] ( PQ{d,H)j^^^ + ^J-^^id,H)Po 



H — En I \ H — El 



^0 



idjH)Po 



H -E, 



^ {d,H)PoPo{d,H) ^ 



H ~ Eq 



1 



H-Eq 



{d,H)PQ 



Tr 



PQ{d,H) 



1 



H-Eq 



{d,H)PQ 



as desired. 



Appendix D: Proof that g is a metric 

By definition, a metric must satisfy three properties |[T]: it must be positive, real, and symmetric. 
(1) Positive: For any nonzero q;(x) e Tm (x) we have 

a(x) • g(x) • q;(x) = gij(x)a*(x)Q;-'(x) 

^-Tr[(9,;Po(x))(a,Po(x))]aXx)aJ"(x) 



250 
Tr 



1 



:a^(x)a,P(x) 



Ik 



1 



a^"(x)a,P(x) 



kl 



= Tr[C^(a,x)C(a,x)] > 0, 



where 



C(a,x) 



:a^(x)a,P(x). 



(C4) 



(C5) 



(Dl) 



(D2) 



Note that although Tr[(dPo)^] is always positive, when we move to a coordinate x the resulting pull-back metric g(x) might 
become singular (non-invertible) at some points or even identically zero. In this strict sense g(x) is not a metric. 

(2) Real: This is obvious from the very construction of g = Re[G]. 

(3) Symmetric: This is obvious from the definition and cyclic invariance of the trace: g^^ = Tr[9iPo9jPo]/2go = 
Tr[9,Po9,Po]/2go = gj«. 



Appendix E: Proof of the operator fidelity inequalities 



We start by proving Eq. (41 1. From the definition of the operator fidelity, Eq. (40 1, with g — I/N, we have, using Eq. ( 14 1: 



Tr 



Tr 



/ 1 



KT{s')n{s')ds' 



> 1- 



nJq 

\Ti[KTn]\ds' 



N 



1 

1 , 

NJo 



Tr([9,,Po,Po]yKd) 



ds', 



(El) 



where in the last line we used the definitions of fi(s) [Eq. ( 13 i] and Kt{s) [Eq. ( 15 i], and cyclic invariance of the trace. Now 
recall the Cauchy-Schwartz inequality for operators 1 18 1: 



IAII2IIPII2 > |(^,i?>| := |Tr[Ats]| 



(E2) 
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Applying this with A := [ds'Po, Pq] and B := VV}^ and noting that |1 VV;'[,|l2 = \/Tr[V;dV"tl/V;|j] = ^fN, we obtain 



f{s)>l-^f \\[ds'Po,Po]hds' = l-^e{s), (E3) 
y N Jo \ N 

as we set out to prove. The inequality /(s) < 1 follows from the fact that Vl{s) is unitary: diagonalizing Vi{s) and taking the 
absolute values of all N of its diagonal elements, which are roots of unity, gives |Tr [fi(s)] | < N . 



Next we prove Eq. ( [44] l. Using Eq. ( 14 1 along with submultiplicativity and the triangle inequality we have: 

\\0{s) - 0.,i{s)\\ = \\0-n{s)^On{s)\\ 

OO CO 

o-{i-Y.^l{smi~j2^^'('y)\\ 



< 



1=1 



l'=l 



OO OO 

1 l'=l 

OO „s 
"'0 



1=1 



o 



1=2 



2 + J2\\n,is 

l'=l 



(E4) 



The term in the first square brackets is identical to that in Eq. (21 1, and hence is bounded by Eq. (22i. The summand 
X^i^i 11^' (^)ll ill f^he second term is 0{1/T) according to Eq. ( 19 1. We thus have: 



||0(5)-0ad(s)|| < \\0\\Ms)+I{s)0{l/T)][2 + 0il/T)] 
where is defined in Eq. ( p5] l, and the last line follows from Eq. ( [T9] l. 



(E5) 



Appendix F: Proof of Eq. j75 



The operator fidelity of two positive operators X and Y relative to a density matrix g is defined as 

feiX,Y)^TT[XYgl 



(Fl) 



which is always nonnegative because the trace of the product of positive operators is nonnegative. When X,Y G GN,ga [subsec 
tion 



IIP 1| and when g is fully supported on the ground eigensubspace, one can conclude from the inequality < Tr[Xy] < 
Tr[y]TT8l that fg{X, Y) < 1. 

Now we compute the fidelity of the ground-state projections -Po(x) and Po(x + dx) relative to g = /go /.go up to the first 
nonvanishing order. Hence, 

/,(Po(x),Po(x + dx)) = (Po(x),Po(x + dx)), 

= -Tr[Po(x)Po(x + dx)] 

= -Tr[Po(x)(Po(x) + dPo(x) + id2Po(x))] 
.90 2 



= 1 



1 

2io 



Tr[Po(x)d2Po(x)Po(x)] 



where in the last two lines we used Eqs. ( Al i and ( A2 1. Equation ( Al i also yields 

d^Po = d^PoPo + 2dPodPo + Pod^Po, 

whence, 

Pod^PoPo = -2Po(dPo)2Po. 



Thus Eq. (F2 1 is simplified as follows: 



1 



/,(Po(x),Po(x + dx)) = 1 - -Tr[Po(dPo)2Po]. 

90 



(F2) 

(F3) 
(F4) 

(F5) 
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Appendix G: Proof of Eqs. {71) and (72) 



To prove Eq. (71 1, we invoke the following inequality 

\Tt[XY]\ < \\XUYl 



(Gl) 



valid for any pair of arbitrary operators X and Y ||57l . In addition, note that by definition, the operator norm of the reduced 
resolvent [H{s) — Eo{s)]^^ satisfies 



His)~Eois) 



< 



dist{{Eoi.s)},spcc{H{.s))\{Eois)}) 
1 



mills ^(-s) ' 

where spec(iJ(s)) is the spectrum of H{s) and the distance between two sets A and B is defined as follows: 

dist(A,B)= inf la -61. 

aeA.beB 



Equation ( 37 1 now yields 



Sij 



< 
< 

Eq.^ 



1 

90 
50 



Tr 



{d,H)Po{d,H) 



1 



H — Eq 



1 



Eq. \G2] , submultiplicativity 



II -Poll 1=30 
< 



1 I 

.go miHs A2 
min„ 



H-Eo 
Po\\i\\d,HdjH\\, 



The proof of Eq. (72i is immediate from noting that |Tr[X]| < J2i (^^W = ll^lli- 



(G2) 



(G3) 



(G4) 



Appendix H: Proof of Eq. i82l 



Note the following identity for the reduced resolvent: 



^OTT ~Qo — Qo 



H — En 



p + H-Eo 



poo 

Qo = / Qoe^-"+"-''°'>^Qo dr 

P=0 Jo 



p=0 



Therefore 



/ 1 \2 d 



1 



dp^p + H-Eo 



Qo = -^ Qoe^-^'+^-^^^^Qodr 
p=o dp Jq 



p=0 



(HI) 



(H2) 



Substituting Eq. (H2i into Eq. (79 1, while recalling that in Eq. (79i the inverse [H — E^] ^ is really shorthand for Qo[H — 
EoV^Qo, yields 



1 d 



90 dp Jo 



p=0 



1 d 
90 dp 7o 
1 d 



dr e-''^Tr[Po{d,Hr)Qoid,H)] 



p=0 



9o dp Jo 



dr e-P^{TT[Poid,Hr)id,H)] - TT[Poid,H)Poid,Hy 



p=0 



(H3) 
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with diHr = e^^diHe-''" . Note that from Eqs. (B8 i and (B13 1, and the property Pq^o^o = 0, we obtain 



Tr[Po(5.ff)Po(a,i/)] - d^EQTr[P^{d,H)] 

= -TT[Po{m)]Tv[Po{d,H)] 
9o 



Substituting Eq. (H4 1 into Eq. (H3 1 and taking the derivative with respect to p yields 

= - / dTTe~P-(TY[Po{d,Hr){d,H)]- -TT[Po{d,H)]TY[Po{d,H)] 
9o Jo ^ 9o 

I 



p=0 



(H4) 



(H5) 



as desired. 



IX. DETAILED DERIVATIONS O F RESU LTS REPORTED 
IN SUBSECTION IIVA3I 

A. Derivation of Eq. jll4^ 

In the thermodynamic limit we replace ^^hy ^^f^ Jq dz, 
where the prefactor is due to a change of variables. Then 
Eq. ( fTTT] l yield 



^ m 



-[l-2(l + C0S5f^)(l-x)x]2 
dz-. 



sin^ z 



[1 - 2(1 + cosz)(l -x)x]2 



2(l-x)2(l-2x)' 
Hence for < x' < 1/2 



(1) 



x(s) 



v/pMdx' = 



x(s) 



dx' 



2 Jo (1 - x')\/l - 2x' 
\ ['^ ^ 4 arctan ^1- 2x(s)] 



1/2 



TT /TT 

2V2- 



(2) 



Now from Eq. ( 1 13 1 we obtain 



2s = 



x(s) 



1/2 



y^dx'/ / VpMdx' 
Jo 

[tt — 4 arctan \/l — 2x(s)] . 
The last equation yields the first case in Eq. ( |114| i, 

x(s) = - l-tan2[-(s- -) 



(3) 



The second case in Eq. ( 1 14 1) is obtained similarly. 



B. Derivation of Eq. jlTT) 

In the thermodynamic limit, using J^e ~^ 
again. 



2m+l 
27r 



, ,.-2cOS:j^l2 
— 1 L 2m+lJ 



dz 



sin^ z 



[1 — 2xcos z + x2]2 



2(1 -x2) 

Hence for < x' < 1 

,x(s) 

V^dx' - 
Now from Eq. ( |1 15| l we obtain 

^x(s) 



— arcsinx(s) 



TT /TT 

2V2- 



yg(xOdx7 / y<z(xOdx' 



arcsinx(s). 



Thus we obtain Eq. ( 1 17 1, 



x(s) = sin(7rs/2). 
C. Derivation of Eqs. (118) and (119) 



To solve 

X + iyKX^/2X = 0, 
(4) (where X = x — Xc) we use the following identity 



x = x^^.lA 

dX 2dX 



X 



/o"dz 



(5) 



(6) 



(7) 



(8) 



(9) 



(10) 
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Hence 



The derivation of Eq. ( 1 19 1 is: 



dX2 

—. = —l^K 



i — 4 X2 = KX-'"' ^ X'^^'/^dX = Kds 



, - = K(s - Sc) 
X{s) = [K{iyK/2 + l)(s - sc)]^"^^ 

= A{.S-Sa)^. (11) 



2 + iyK 

Eg . (87) 2 

~ 2 + v{2a -2z-d) 

a—d+z—l/u 2 



Therefore 



2 + v{2d +2z- 21 V -2z-d) 
2 

dv' 



(13) 



x(s)-Xc = A(s-Sc)2+-'=. 



(12) 
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